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Abstract 



. We study the chiral phase transition in the hnear sigma model with 2 quark flavors 

' and Nc colors. One-loop calculations predict a first-order phase transition at both 

lO ' = and ^ 7^ 0. We also discuss the phase diagram and make a comparison with a 

thermal parametrization of existing heavy-ion experimental data. 

^ : 1 Introduction 

p • The linear sigma model, originally proposed as a model for strong nuclear interactions 0, 
^ ■ today serves as an effective model for the low-energy (low-temperature) phase of quantum 
chromodynamics. The model exhibits spontaneous breaking of chiral symmetry and restora- 
tion at finite temperature. In this paper we discuss the pattern of symmetry breaking and its 
^ I restoration at nonzero temperature and nonzero chemical potential at one-loop level. Some 
c3 ' aspects of this have been extensively studied in the literature [2-12]. The phenomenological 
importance of the chiral transition and possible experimental signatures have recently been 
discussed by Harris and Miiller . In connection with theoretical predictions of the phase 



transition there still exist a certain misunderstanding and controversy. In particular, the pre- 
cise value of the critical temperature is not generally agreed on [0, |^, |^, . As regards the 
order of the transition, Pisarski and Wilczek have shown on general grounds that chiral 



models with 2 massless flavors undergo a first-order phase transition at finite temperature. 
In actual calculations, the mean field predicts a second-order transition P, whereas some 



improved mean-field results indicate a weak first-order transition Numerical simulation 
||T5| of a three-dimensional chiral model on the lattice also confirms a first-order transition 
for Nf = 2. We believe that in our approach we are able to clear up some of these points. 
In addition, we discuss the phase diagram and its relation to heavy-ion collisions. 

We organize the paper as follows. In Sect. |^ we describe the thermodynamics of the 
model. In Sect. we calculate the temperature dependence of the pion and sigma masses. 
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In Sect. ^ we discuss the chiral phase transition and make comments on and a comparison 
with other papers. In the concluding section, Sect. ^, we summarize our results. 



2 Effective potential 

The linear sigma model of Gell-Mann and Levy is an effective model of strong interactions 
described by the chirally symmetric Lagrangian (in Euclidean notation) 

£ = tp^j^d^ + g{a + iT-n-'j^))^ - fitpjii; 

+ l{dar + lidT^r + ^{a^ + TT^) + ^{a^ + tt^)^ . (1) 

The right and left fermions ipR = 1(1 + 75)"^) i'L = |(1 — 15)4', are assumed to constitute, 
respectively, the (^,0) and (0, |) representation of the chiral SU(2) xSU(2), whereas the 
mesons (a, tt) belong to the (|, |) representation. In the original sigma model [|l| the 
fermion field was a nucleon. We consider the fermions to be constituent quarks ^ with an 
additional degree of freedom, "color" , from the SU(Ai'c) local gauge group of an underlying 
gauge theory (QCD). 

If rriQ < 0, the chiral symmetry will be spontaneously broken. At the classical level, the 
a and vr fields develop nonvanishing expectation values such that 

{ar + {7rr = -^^f.. (2) 

It is convenient to choose here 



(tt,) = 0, (a) = U ■ (3) 

In order to study the thermodynamics of the model, we define the thermodynamical 
potential as a function of the chemical potential /i associated with baryon number density 
and inverse temperature f3 = 1/T 

n{(3,^^) = -^\n Z , (4) 

where the partition function Z is defined as a path integral 

Z = J [dip] exp{- J^^ dr J d^x C{ifi)} . (5) 

Here [dip] is an abbreviation for the integral over tp, a, tt, and C is given by (|l|). Next we 
introduce the saddle-point method of Frei and Patkos [T^ . Our approach is similar to that 



of Meyer-Ortmanns and Schaefer [|T7| who applied the method to the chiral SU(3)xSU(3.) 
We first redefine the fields 



TT — i> 77 + 7r'(x) , 



a ^ a + a'{x) , (6) 
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where tt' and a' are quantum fluctuations around the constant values tt and a, respectively. 
Next we use the transformation which quadratizes the quartic interaction 

/\ /-c+ioo fa 

d'x \{a'^ + Tr'^)^ = \ds\ exp{- / d'x ^-{a'^ + tt'^) - (7) 

and redefine the auxiliary field 

s{x) ^ s + s\x) , (8) 

so that the saddle-point value s maximizes the integrand. The thermodynamical potential 
as a function of cr, tt and s is usually called the effective potential. Thermodynamics requires 
that, in thermal equilibrium, f2 should assume a minimum with respect to variations of tt 
and a . Owing to the specific form of the interaction and because of (|^) and (^), we can keep 
77 = and consider as a function of the two mean fields o and s. Using (^, (|^) and (|) 
we find 

^{o. ^) = I + ^ - Ix - ^ ' 

where Z' is the partition function for the shifted Lagrangian in which the quartic interaction 
is absent and chiral symmetry is explicitly broken: 

a = ^(7(9 + mF + c/(a' + iT7r'75))^-/xV^74V^ + ^((97r)'2 + l(5(T')^ 

+ ^a'' + ^TT- - + gW{a- + tt'^) + '-{a'' + tt^) + ca' - l-ss' . (10) 

The effective masses, the trilinear coupling g' and c are functions of a and s defined as 

m\ = ttIq + s + 3A(j^ , nip = ga , 



rni = mf. + s + Acr^ , q' = \a . 



c = a{ml + Xa^) . (11) 
The condition for an extremum 

gives two equations for a and s with solutions that will in general depend on temperature 
and chemical potential. A nontrivial solution a{f3, /z) will be referred to as chiral condensate. 

At the classical level (neglecting the quantum and thermal fluctuations) the potential 
takes the form 

Q{a,s) = lxa' + ^a'-^. (13) 
An extremum (minimum with respect to a, maximum with respect to s) occurs at 

A 

yielding 



2 



a 



-^^f., s = 0, (14) 



0,mF = g^,ml = 2\fl , g' = Xf^ , c' = . (15) 
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Thus, at the classical level, we have 3 massless pions, as it should be owing to the Goldstone 
theorem. 

The thermal and quantum fluctuations will change the effective potential into 



(16) 



Here Qj contains loop corrections and Qq is the thermodynamical potential for a noninter- 
acting gas of fermions and bosons: 



fig — -|- fig- + 



(17) 



E 



2P^J (27r)3 



In 



where 



TC 



p= ( (2/+ l)- + z/i,p 



TT 



k = ^n- , kj , 

and the number of flavors is Nf = 2. The extremum condition now reads 



(18) 

(19) 
(20) 



on 



1 



Xa^ + mla + j d^x { g{tlj{x)^{x)) + 3Xa{a'{xy) 



+\a{7T'{xY) + \{a\x){a'{xY + 7z\xY)) + {ml + 3\a'){a'ix)) = , (21) 



(22) 



The terms {a') and {s') vanish because the quantum fluctuations take place around the true 



vacuum. Equations (0) and (^2]) become 

Xa^ + mla - 2gN,gF + 3Xa[g^ + G^] + A[r,,, + T^^^] = , 

Here Q and F denote thermal averages (with respect to the full partition function) over a 
product of two and three flelds, respectively. For example. 



(23) 
(24) 



= ^ / d'x {a{xf) . 



(25) 



This coincides with the full 2-point Green's function at a; = 0. ^ is often referred to as 
tadpole. Equations (|23| ) and (|2^ , schematically depicted in Fig. 0, state the fact that 
tadpoles cancel |]18| also at nonzero temperature and chemical potential. 

Solutions to (p3D and (|^) are implicit functions of T and fi. The tadpoles at one-loop 
order are given by 



Qf 



Tr 



1 



/3 



E 



—if + mp 
1 



(27r)3F + m2^ ' 



(26) 



whereas the three-point functions Fo-o-o- and Fo-tttt contribute at two-loop (and higher) order. 
The masses in (p6D depend on a and s through ([TT|). Equation ( ^3]) has, apart from cr = 0, 
a nontrivial solution a{P,fi) that no longer equals /,r- Dividing it by Act, this equation at 
one-loop order may be written as 



d^k 



(27r)3 p'^ + rup 



3iE 



(27r)3 F + m2 (3 



d^k 



(27r)3p + m2 



where we have replaced — mg/A = from (0) and used mp = go-- Similarly, 



1 r d^k 1 1 r 

/3 ^ 7 (27r)3 P + m2 ^ /3 ^ 7 



:27r)3fc2 + m2 



(27) 
becomes 
(28) 



We can separate the finite T and /i part of Q as usual |19 

d^q 



mp 



and 



where 



(27r)4 g2 + ^2 
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(27r)4 g2 + ^2^ 



1 

ni?(a;) 



+ 



(/■^g mp 
(27r)3 2Z7^ 

d^g 1 



2 I 2 
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(29) 

(30) 

(31) 
(32) 

(33) 

The infinite part in (pQf ) and ( pOD is the usual T = ^ = tadpole that is absorbed in the 
tadpole cancellation at the tree level. Therefore, we can write our equations by retaining the 
T- and /^-dependent pieces only. Equations (p7|) and (pHf) finally read 



+ 



1 
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a 



f d^q 1 f d^q 1 



A 



(27r)'^ o'o- J (27 

where the right-hand sides depend on a and s through the masses. These equations have 
been derived from the effective potential (|16]) in which the loop corrections VLj have been 
neglected. This approximation corresponds to the leading order in the expansion, where 
N is the number of scalar fields Jl^. In our case, = 4. 

A straightforward approach to solving ( |34D and (^) leads to problems with a complex 
effective potential [jl^. It may easily be seen that a direct use of (|ll]) leads necessarily to 



complex solutions. From (0), (^) with ([TT]) one finds 

2 2 f d^q 1 . f d^q 1 , , 

= -8, AT, y ^ ^ - 2A / ^ - usM- (36) 

This means that is either negative or complex. In both cases it implies a complex a. In 
the following sections we show how a consistent inclusion of one-loop self- energy corrections 
removes this problem. 

3 Effective meson masses 

At first sight it seems that the T, /i dependence of a and s is not consistent with the Goldstone 
theorem since 

ml + s{/3,n) + Xa\f3,fi)^0 . (37) 

However, the must also include the T- and fi- dependent pieces coming from the one-loop 
(and higher) order self-energy diagrams (Fig. ^): 

ml = ml + s{(3, /i) + Xa\l3, /x) + n^(/3, /x) , (38) 

where (at one-loop order) 



n = -8Nq'-Tf^^ ^q"-Tf 



2 



(39) 

with p and k defined as in (|l^) and (|20|) , respectively. The s-field appears in the Lagrangian 
without a kinetic term and its propagator is simply l/m^ = —2 A. 
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The combined contribution of the a-tadpoles (pTj) and the s-tadpoles ( p8D is given by 




(40) 



Combining (|39D and (0), we see that the fermion parts cancel immediately, whereas the 
boson parts cancel provided the relation 



which holds trivially at T = = 0, also holds at T, < Tc, Hc- Thus, the consistency with 
the Goldstone theorem requires that the relation 



should also hold at T, /x < Tc, /ic- Indeed, we shall shortly demonstrate that it works at 
one-loop level. 

We do not agree here with Larsen who obtained ^ even in the symmetry-broken 
phase, thus violating the Goldstone theorem. The reason is that he assumed m^^ 7^ in the 
propagators in the one-loop self-energy diagrams. If = at the tree level, then the mass 
corrections would be of order A or higher. This would in turn yield self-energy corrections of 
order A^, which would then also require the inclusion of two- loop diagrams. In other words, 
it is not consistent with one-loop calculations to include corrections to the tree-level masses 
in the propagator. 

Similarly as for the pion mass, for the sigma mass we have 




ml = 2Xa^ , 



(41) 



ml = 2\(j^ 



(42) 



ml = ml + s{P, /i) + 3Acr^(/3, /i) + n^(/?, /i). 



(43) 



where the self-energy of the a particle (Fig. |^) is given by 




d^k 1 

(27r)3 (P + ml) 



d^k 1 

(27r)3 {k^ + miy 




(44) 



This, together with the tadpole part 




2A/2 




(45) 
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gives 



ml 



P^J (27r)3 p^ + ml 



I f d'^k 1 I r d'^k 1 

^ ^^^^^J (2^ k^ + ml~ ^ (2^ (F + m2) ' ^^^^ 

where ... denotes terms of higher order in A and g^. Combining this equation with (|27|) , 
we recover (^). Thus, the thermal corrections at one- loop level do not alter the tree- level 
equation for the effective sigma mass. 



4 Chiral restoration transition 



From the analysis in the preceding section we see that the s-dependence of the rrifj and m^r 
is removed owing to the one-loop self-energy corrections. Equations (|5^, (^) in which we 
put 



ml 
mp 



0, 

, 



(47) 



decouple and we are able to determine the T and n dependence of the chiral condensate by 
solving (0) only. If one compares our approach with the standard one where the quartic 
interaction is kept, one may easily check that the self-energy diagrams with internal s-lines 
plus the contribution of s-tadpoles precisely equals the contribution of the loops with the 
quartic vertex. In this way, the saddle-point method [T^, ^ becomes equivalent to the 
standard approach. 

For /i = 0, we find the solution numerically as a function of T depicted in Fig. |^. As 
input parameters we choose the constituent quark mass mp = 340 MeV and the sigma mass 
m^j = 1 GeV. The solution indicates a first-order phase transition. The point where the 
curve crosses the T axis is not the point of the actual phase transition. We refer to it as 
"critical" point, having in mind that it would be a critical point if the phase transition were 
second order. The actual transition takes place at the point where the two minima of the 
effective potential at a = and ct(Tc) are leveled. Hence, the transition temperature Tc is 
determined by requiring 

n{a,T,) = n{0,T,) . (48) 

The solution is particularly simple in the neighborhood of the "critical" point since in 
this case we can expand the integrands around a = and perform the integrals analytically 
g. We find 



J n 



— + — 
A 6 2 



+ 



3V2A 



47r 




2g^N, 
A7r2 



+ 



3A 



In 



2Aa 



47r2 47rT 



g'Nc^^ga 



+ 0{a^), (49) 
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where 7 = 0.5772... is Euler's constant. At the "critical" temperature, cr(T^') = , which 
gives 

T'2 6A 



/2 2^?2iV, + 3A 

If we had taken the Nambu-Jona-Lasigno relation ma = 2mF [pOl' which is exact in the 



(50) 



iVc — >■ cxD limit we would have obtained T^' = /j^. We stress again that only approxi- 
mates the actual transition temperature to be determined from the condition 

If we put Nc = 1 vci (|50D our result agrees with that of Anand et al. 0] who considered 
a similar model for nuclear matter in Walecka's mean-field approach. Our result (^) also 
agrees with Bochkarev and Kapusta [|12| if we put = 0. However, some of the existing 
calculations that include fermions 0, |^, ^ disagree with our result for the following reasons. 
Ram Mohan p[, and subsequently Contreras and Loewe did not properly account for the 
contribution of antifermions. As a consequence, they obtained fewer (by a factor of two, 
apart from color) fermionic degrees of freedom, which in turn yielded a larger estimate for 
the Tc The reason for disagreement with Cleymans, Kocic and Scadron is twofold. 
First, in their calculation of the sigma mass they did not include all the relevant self-energy 
diagrams (Fig. Second, the sign of the fermionic contribution to the sigma mass is 
wrong. This led again to a larger in terms of /t^. The agreement of their estimate of the 



critical temperature with the estimates based on a single meson- loop diagram [11] and on 



the Nambu-Jona-Lasigno model is only accidental. 

The calculations at nonzero chemical potential are similar. We present our results in 
Fig. The chiral condensate a is plotted as a function of temperature for fixed fi (160 
and 350 MeV; upper two plots) and as a function of chemical potential for fixed T (0 and 
50 MeV; lower two plots). The transition remains first order and the critical temperature 
decreases with fi as expected. 

In Fig. P we plot the phase diagram of nuclear matter and compare it with the thermal 
parametrization of recent heavy-ion collision data. The baryonic chemical potential is related 
to the quark chemical potential as /i^ = 3fi. The phase boundary between the chirally 
symmetric and broken phases (solid line) appears to be very close to the expected phase 
boundary between hadron resonance gas and quark-gluon plasma with a bag constant [|13|. 
The diagram indicates that the chiral transition might slightly precede the deconfinement. 
The points with error bars show the freeze-out values of T and /i deduced from AGS [^, ^ 
m and SPS ^ data with flow. 



In Fig. ^ the phase diagram is represented in terms of temperature and baryon density. 
The solid line separating the chirally symmetric and broken phases shows how the transition 
temperature depends on baryon density. 

Similarly to the chiral condensate, the effective masses will have a discontinuity at Tc. 
Above the critical line in Fig. ^ the symmetry is restored and the a-tadpoles vanish from 
the theory. The T and fi dependence of the meson mass itlm = = is determined by 
the proper self-energy diagrams and s-tadpoles: 

m^^ = mo + s(/5,/i) + n^(/5,/i) . (51) 



Here s{P,fi) and may be calculated at lowest order in A and using (pSD and (^ 

in which the propagator masses and g' are set to zero. In particular, for /i = and above 
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we find the following expression for the meson mass 

-l.= (^ + f/)(T^-Tf). (52) 

5 Conclusion 

We have shown that the usual mean-field pattern of symmetry breaking and restoration 
gives a consistent picture in the cr-model at the one-loop order. We have shown that the 
saddle-point method [jl6|, |l^ is equivalent to a standard approach ||, ^, |T2| if the self- 



energy loop corrections are included. The phase transition is predicted to be first order, 
in agreement with the analysis of Pisarski and Wilczek [Q. The chiral phase boundary 
in a (T, /is) plot (Fig. ^ is close to the phase boundary between hadron gas and quark- 
gluon plasma. The thermal parametrization of existing experimental data compared with 
the chiral phase diagram indicates that the nuclear matter produced in heavy-ion collisions 
is close to or slightly above the chiral-phase-transition line. It is therefore conceivable that 
the present and future heavy-ion experiments may observe effects of the chiral transition. 
A more sophisticated analysis of the data is needed in order to observe possible signatures 
near the critical density, such as a rapid change of the meson mass and width or abnormal 



production ratios of charged to neutral pions ||13 . 
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Figure 1: Schematic representation of (23) and (24) 



Figure 2: One- loop self-energy diagrams contributing to (a) pion and (b) sigma masses. 



Figure 3: Chiral condensate as a function of temperature at // = 0. The dashed line cor- 
responds to a physically unstable solution. is the temperature of the first-order phase 
transition. 



Figure 4: Chiral condensate at nonzero baryon density, ex, T and ji are in MeV. 



Figure 5: Phase diagram as a plot of temperature versus baryon chemical potential. The 
solid line separates the chirally broken (inside) and chirally symmetric (outside) phases. The 
expected hadron gas - quark-gluon plasma boundary is located between the two dashed lines. 



Figure 6: Phase diagram as a plot of temperature versus baryon density. 
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